Hydrodynamic coupling of two rotating spheres trapped in harmonic potentials 
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We theoretically study in detail the hydrodynamic coupling of two equal-sized colloidal spheres at low 
Reynolds numbers assuming the particles to be harmonically trapped with respect to both their positions and 
orientations. By taking into account the rotational motion, we obtain a rich spectrum of collective eigenmodes 
whose properties we determine on the basis of pure symmetry arguments. Extending recent investigations on 
translational correlations [J.-C. Meiners and S. R. Quake, Phys. Rev. Lett. 82, 221 1 (1999)], we derive the com- 
plete set of auto- and cross-correlation functions emphasizing the coupling of rotation to translation which we 
illustrate in a few examples. An important feature of our system is the self-coupling of translation and rotation 
of one particle mediated by the neighboring particle that is clearly visible in the appropriate auto-correlation 
function. This coupling is a higher-order effect and therefore not included in the widely used Rotne-Prager 
approximation for the hydrodynamic mobilities. 

PACS numbers: 82.70.Dd, 05.40Jc, 47.15.Gf 
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I. INTRODUCTION 

Colloids are widely used to model atomic systems (H • 
However, there is one feature specific to colloidal suspen- 
sions which distinguishes them fundamentally from atomic 
systems: the so-called hydrodynamic interactions |3l|4[. A 
particle moving in a viscous fluid creates a long-range flow 
field around itself through which it interacts with other par- 
ticles. Thus, hydrodynamic interactions constitute a compli- 
cated many -body problem since the motion of one particle de- 
pends on the translations and rotations of all the other particles 
in the fluid H. 

The central quantities in the description of hydrodynamic 
interactions are the mobility or friction tensors which connect 
in a linear response scheme all the forces and torques acting 
on the particles to their linear and angular velocities [6]. In the 
present paper, we draw special attention to the rotational de- 
gree of freedom and how it couples to translation. Itsphysical 
consequences have rarely been treated in literature QIMIal- 

In many physically interesting systems, hydrodynamic in- 
teractions play an important role. Whenever dynamic effects 
in suspensions are studied, hydrodynamic interactions have 
to be taken into account. Examples are sedimenting particles 
ifTot Hill or the apparent attractive interaction between like- 
charged spherical particles mediated by a single, like-charged 
wallTTl. 

In recent years, a new method called microrheology has 
been dev elop ed I13ll and applied to several biological sys- 
tems lfl4l fl5fl . It is used to determine rheological properties 
of viscous and viscoelastic media by tracking the trajectory of 
embedded probe particles and calculating the time-dependent 
position correlations. The so-called two-point microrheology 
employs a system of two colloidal particles, which has sev- 
eral advantages over using only a single particle; in particu- 
lar, the correlated motions of two tracer particles reflect the 
bulk rheology of the medium they are embedded in more ac- 



curately fl6lfl7ll . Many experiments based on this technique 
have been carried out to study soft media fT8lfT^.l20ll2lll . Fur- 
thermore, this method was used to measure the corrections to 
the diffusion coefficients in a system of two spheres due to hy- 
drodynamic interactions 1 22] and to study the hydrodynamic 
coupling of two spheres to a wall 

The physical systems we have in mind are colloidal suspen- 
sions of birefringent spherical particles. E.g., they are realized 
by p olymerizing the nematic order in liquid-crystal droplets 
11241 12511 . In contrast to conventional isotropic colloids, both 
their positions and orientations can be manipulated by opti- 
cal tweezers. In particular, the orientation is controlled with 
optical traps generated by linearly polarized laser light and is 
detected, e.g., with crossed polarizers (see Fig.^. 

In the present paper, we consider a system of two trapped 
spheres of equal size. We give a complete analytic descrip- 
tion of the coupled translational and rotational dynamics. Due 
to the axial symmetry, the longitudinal motions parallel to 
the particle-particle axis decouple completely from transver- 
sal motions. In addition, there is no coupling of translation 
and rotation in longitudinal motions. This greatly reduces the 
complexity of the problem. We determine the complete set of 
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FIG. 1 : Sketch of a setup to monitor rotational correlations of two 
birefringent particles. The particles are situated on the z-axis. The 
laser beam of the orientational optical trap is directed along the y- 
axis and its polarization points into x-direction. To monitor transver- 
sal correlations, the direction of observation is along the y-axis. Pure 
longitudinal rotational correlations are only observable along the z- 
axis. 
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collective eigenmodes of the system and calculate the correla- 
tions for particles undergoing thermal fluctuations. Our work 
extends recent investigations on correlations in translational 
fluctuations lITsL IT9II2OII . By taking into account the rotational 
degrees of freedom, we obtain a rich spectrum of collective 
modes and correlation functions. 

The outline of this paper is as follows. In Sec. [H] we first 
summarize the basic equations for the description of hydro- 
dynamic interactions in general. Then, we study in detail the 
two-particle system and show by symmetry arguments how it 
decouples into the aforementioned longitudinal and transver- 
sal motions. Assuming that both the particles' positions and 
orientations are trapped in harmonic potentials, we determine 
and discuss their collective eigenmodes in Sec. [TO] Finally, 
in Sec. lIVI we use Langevin dynamics to calculate correlation 
functions for the thermally fluctuating positions and orienta- 
tions of the trapped particles and point out their interesting 
features. 



n. HYDRODYNAMIC INTERACTIONS 

A. Fundamental equations 

In the regime of low Reynolds numbers and on the Brow- 
nian time scale, the flow of an incompressible fluid with vis- 
cosity 77 obeys the Stokes or creeping flow equations |3], 



r)V 2 u - Vp = , V 



0, 



(1) 



where u is the flow field and p the hydrodynamic pressure. 
Stokesian dynamics describes overdamped motion in a vis- 
cous fluid. In the following, we consider motions of particles 
in an unbounded and otherwise quiescent fluid, i.e., u = at 
infinity. 

Imposing stick boundary conditions on the surfaces of all N 
particles suspended in the fluid at positions r» (i = 1 . . . , N), 
the motions of the particles are mutually coupled via the flow 
field. Due to the linearity of Eqs. 0, the translational and ro- 
tational velocities of the particles, V{ and u>j, depend linearly 
on all external forces and torques acting on the particles, Fj 
andTj H: 



N 



0) > 



N 



i=i 



(2a) 
(2b) 



The central quantities constituting the mutual coupling of 
translation and rotation (denoted by superscripts t and r) of 
two particles i and j are the 3x3 mobility tensors /i,**, fxlj, 
fj,\j, and fi 1 ^ . They depend on the current spatial configuration 
of all particles, i.e., the set of position vectors {t*i, ...,»•#} 
in the case of spherical particles. 

To introduce a more compact notation, we define the 6N- 
dimensional vectors V = [v\, . . . ,vn,lji, . . . ,u>n} andF = 



[Fi, . . . , F N: Ti, 



, TV]- Then, Eqs. (0 take the form 
V = MF (3) 



with the 6iV x 67V mobility matrix 



M 



(4) 



where the four blocks of M consist each of N x N matrices 
whose elements are the 3x3 mobility tensors; e.g., 



tt 

A* 11 



tt 
r*lJV 



t^Nl ■ ■ ■ t^NN 



(5) 



According to the reciprocal theorem of Lorentz, the mobility 
tensors fulfill the symmetry relations |6] 



(M") T 



(6) 



where (/^") T denotes the transpose of /iy, etc. Thus, the 
entire 67V x 6iV matrix M is symmetric. 

In the overdamped limit discussed here, the work done on 
the particles by the external forces and torques F is com- 
pletely dissipated in the fluid, so the energy dissipation rate 
■ Fi + u>i ■ Ti) = V • F = V • M^V has to be 
positive. Therefore, the friction matrix M _1 and hence the 
mobility matrix M itself are positive definite. 

There a several methods to calculate the mobility tensors 
for a given many-particle system, e.g., the concept of reflected 
flow fields using the gradient expansion technique 1 3 , 6 ] or the 
method of induced force multipoles The latter was 

implemented in the numerical library HYDROLIB |28] which 
calculates the mobility or friction matrix for a given configu- 
ration of spheres. 



B. Two-sphere system 

In the following, we consider a system of two equal-sized 
spheres. Let r be the vector connecting the centers of the two 
spheres, pointing from sphere 1 to sphere 2, and r the center- 
to-center distance. Then, r = r/r is the unit vector along the 
line of centers. Due to the rotational symmetry about the axis 
r and the different parities of polar (vi and F j) and pseudo- 
vectors {u>i and Ti), the mobility tensors can be written as 
1291 



(4j (r) r <g> r 



f ® f ) , 
r <g> f) , 



(r)= M *f(r)fx 



(7a) 
(7b) 
(7c) 



The mobility coefficients l^j'^, an d lAj 1 ' (h 3 = 1; 2) 

are scalar functions depending only on the center-to-center 
distance r. They describe motions parallel and perpendicular 
to the axis, respectively. Using the general symmetry relations 
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FIG. 2: Definition of the applied Cartesian coordinate system. The 
z-direction is along the center-to-center line of the two spheres. Lon- 
gitudinal motions (a) decouple from transversal motions (b). Further- 
more, there is no coupling of translation and rotation for longitudinal 
motions. 



and the fact that particles 1 and 2 are identical, the mobility 
coefficients obey [29] 



ttll,_L _ ttll,_L 
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(8a) 
(8b) 



With the explicit expressions Q, the dynamics of the two- 
sphere system described by Eq. separates into longitudinal 
motions parallel and transversal motions perpendicular to the 
particle-particle axis f. Furthermore, longitudinal translations 
and longitudinal rotations (i.e., rotations about the axis) do 
not couple to each other because of the different parities of 
translations (polar vectors) and rotations (axial vectors). In 
a Cartesian coordinate system with the z-direction pointing 
along f (see Fig. [2), we obtain from Eqs. 10, Q, and (|8} for 
longitudinal motions 



(9) 



where v\ z is the z-component of the translational velocity of 
particle 1, etc. In transversal motions, translations along the 
x-direction are coupled to rotations about the y-axis: 
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(10) 



Translations along the y-direction and rotations about the x- 
axis obey an equivalent system of equations. 

Defining a four-dimensional velocity vector v and a force 
vector f, we abbreviate Eqs. (|9jl and i ll Oi l, respectively, by 



mf . 



(11) 



where the appropriate 4x4 mobility matrix m is still sym- 
metric. Thus, rotational symmetry reduces the full 12 x 12 
problem Q to two 4x4 problems, where Eq. (|9jl is essen- 
tially a 2 x 2 problem. 

Accordingly, the energy dissipation rate V ■ M _1 V splits 
up into a sum of three independent terms of the form v ■ m 1 v, 
so the single mobility matrices m have to be positive definite. 



C. Mobilities in Rotne-Prager approximation 

For a system of two particles, the mobility tensors can be 
calculated, e.g., with the method of reflections based on the 
Faxen theorem. It provides a systematic expansion of the mo- 
bility tensors in powers of the inverse particle separation (for 
a detailed description, see, e.g., Ref. y|])- 

The leading order in the far-field approximation of the mo- 
bilities is the well-known Oseen tensor which considers the 
particles as point-like and hence does not include rotations. 
The next-higher order is the so-called Rotne-Prager approxi- 
mation. It corresponds to one reflection of the flow field and 
is therefore exact up to the order of 1/p 3 , where p = r/a 
is the dimensionless particle distance (a is the radius of the 
spheres). The relevant mobilities including the rotational de- 
grees of freedom are presented, e.g., in Refs. l29l l3(ill . They 
also follow straightforwardly by extending the calculations in 
Ref. 1 3], where the torques are explicitly set to zero. 

In the Rotne-Prager approximation, the self-mobilities are 
identical to the Stokes coefficients for single spheres, 



M, 

m! : 



m'i- 

M T 1: 



where 
where // 



(12) 
(13) 



and there is no self-coupling of translation and rotation, i.e., 



t4 = t4 = o. (i4) 



Note that for an isolated sphere within the linear Stokes 
regime, translation and rotation are not coupled; the so- 
called Magnus effect only enters via the nonlinear term in the 
Navier-Stokes equation 1 3 1 ] . On the other hand, in the two- 
sphere system, the self-coupling exists since it is mediated by 
the second particle. However, this is an effect of higher or- 
der than 1 / p 3 and therefore not included in the Rotne-Prager 
approximation, as stated in Eq. dl4l >. 

For completeness, we also give the essential cross- 
mobilities: 



M12 — M 



rr r 
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, (15) 
(16) 
(17) 



III. EIGENMODES OF TWO TRAPPED SPHERES 

We aim to study thermal motions of particles that are 
trapped with respect to both their positions and orienta- 
tions. Assuming that the spatial and angular displacements 
are small, we consider harmonic trap forces and torques, i.e., 



F l( 
F„ 



ICt 1 

ia. j 



(18) 
(19) 
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FIG. 3: While the rotation of the optical axis of the trapped particle 
about the direction of the trap polarization (here pointing along the 
x-axis) is free, the deviation from the ^-direction is restricted, and 
therefore the angular displacements \iy an d X»z are small (see text). 



where fc* and k r are positive force and torque constants (i = 
1, 2 is the particle number and a = x, y, z the coordinate in- 
dex). The particle velocities are 

Via = ha , (20) 

where the spatial coordinate Ti a denotes the displacement of 
particle i along direction a. 

The angle x%a describes the rotation of particle i about the 
coordinate axis a. As long as this angular displacement is 
small {(Xia) ^ 1)' we can relate it to the angular velocity via 



Xi 



(21) 



This is a nontrivial statement and requires some explanation. 
The optical axis of the trapped birefringent particle aligns 
along the polarization of the laser trap pointing along, e.g., 
the x-axis (see Fig. [3}- Any deviation of the particle ori- 
entation from the x-direction is described by the angles Xiy 
and Xiz, and relaxes in the optical trap on the time scale 
(fe r //) _1 . However, the free rotational diffusion of the par- 
ticle axis about the x-direction also changes the angular dis- 
placements Xiy an d Xiz (see Fig. [3}, and therefore, relation 
i21\ does not hold in general. Nevertheless, in the limit con- 
sidered here the orientational relaxation is a much faster pro- 
cess than the free diffusion, so on the relaxation time scale 
(fc r //) -1 , the free diffusion can be neglected and relation i2l\ 
is applicable. To show this, we note that the free rotational 
diffusion takes place on the time scale (feeT^ 1 ) -1 . Requir- 
ing that OV^VObTX)- 1 = ksT/k 1 < 1 and using the 
equipartition theorem yields (xf Q ) = k-^T/k 1 ' <C 1, which is 
consistent with the initial assumption. 

We define generalized coordinates q = [n x , r2 Z , Xiz, Xiz\ 
for longitudinal modes and q = [r\ x , r2 X ,Xiy> X2y] for trans- 
versal modes and combine the corresponding trap forces Jl 8i 



and torques dl9t into a four-dimensional vector 
f = kq 

with the diagonal force constant matrix 



fc* 

fc* 

A: 1 ' 

k 1 



(22) 



(23) 



Then, the equation of motion (^J for both longitudinal and 
transversal modes reads 



q + mkq = . 



(24) 



where m is the appropriate mobility matrix given by Eqs. ij^Ji 
and il OK respectively. Note that for spatial displacements 
small compared to the equilibrium particle distance, we can 
consider m to be constant. 

The solutions of Eq. 1241 are relaxational eigenmodes 
e "*a„ with relaxation rates X n or relaxation times A^ 1 . 
They are determined by the eigenvalue problem 



mka„ = A„a„ (n = 1, . . . , 4) 



(25) 



of the nonsymmetric 4x4 matrix mk whose eigenvectors a„ 
are in general not perpendicular to each other. In the next two 
sections, we will analyze the eigenmodes in detail. 



A. Longitudinal eigenmodes 

To study the longitudinal eigenmodes, the traps have to be 
polarized along the x- or y-direction so that the angular dis- 
placements Xiz are small. Since longitudinal translations and 
rotations are decoupled (as discussed in Sec. Ill B> . we can im- 
mediately write down the eigenvalues and eigenvectors of the 
four eigenmodes: 

A1/2 - ** (Mn ± Ma') , a 1/2 = [1, ±1, 0, 0] , (26a) 
A3/ 4 = fe r (Mi r i"±Mi2% a 3 / 4 =[0,0,l,±l]. (26b) 

They consist of relative (— ) and collective (+) modes. In- 
tuitive arguments allow a comparison of the respective re- 
laxation rates. E.g., when the spheres translate in opposite 
directions (relative translational mode a 2 ), some fluid has to 
be pulled into or squeezed out of the region between the two 
particles; or when they rotate in opposite directions (relative 
rotational mode aj), the fluid between the spheres has to be 
sheared. On the other hand, when the spheres translate or ro- 
tate collectively (modes ai or a3), the fluid surrounding the 
spheres is just "displaced" or "rotated" as a whole. So the 
collective modes experience less resistance and therefore re- 
lax faster than the relative modes. This is in accordance with 
Eqs. J26l > which give Ai > A2 and A3 > A4. 
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B. Transversal eigenmodes 

The appropriate coordinates to treat the transversal eigen- 
modes are q = [r Xx , r 2x , Xiy, Xiy\, and the mobility matrix 
m is given by Eq. dlOi . The polarization of the traps has to be 
chosen along the x- or z-direction so that the angles \iy are 
small. 

The diagonalization of the nonsymmetric 4x4 matrix mk 
is not immediately obvious. However, symmetry arguments 
help to identify the eigenvectors. Since the two particles are 
identical and since spatial coordinates and angles possess dif- 
ferent parities, one readily shows that the eigenvectors have to 
be of the form [A, A, B, —B] or [A, -A, B, B]. This con- 
straint simplifies the determination of the eigenmodes con- 
siderably. There are two relaxational modes with symmetric 
translation and antisymmetric rotation, 



A1/2 



fc>+ + fc> (27a) 



a l/2 = [A-xl2,A x i2,B x /2,—Bx/2 



1/2 

where 

^1/2 
B X /2 



(27b) 



-2)fc> t _!- L , (27c) 

-(feV? X -fc r M- X ) (27d) 

± \J (fc'/i"- 1 - fc 1 '/^) 2 + 4fc t fc r (/x t _ r - L ) 2 , 

and two modes with antisymmetric translation and symmetric 
rotation, 



A 



3/4 



a 3/4 = [^3/4, -^3/4, -B3/4, -B3/4] 

where 



.4 



3/4 



B 



3/4 



(28a) 



(28b) 

(28c) 
(28d) 



± J (fcV U± - fc 1 '/^) 2 + 4fc t fc r (^ a ) 2 



Here, we use the abbreviations /z^ = ± 1 > etc. 

Without relying on the exact values of all the mobilities, we 
can already infer the qualitative characteristics of the eigen- 
modes based on general arguments. The signs of the compo- 
nents Bx-4 are independent of the mobility coefficients. It 
is Bx > 0, B 2 < 0, B 3 > 0, and B A < 0. The signs of 
Ax-a depend on the signs of /J t ± ± - In Sec. Ill CI we already 
discussed on the basis of the reflection method that the abso- 
lute value of the cross-mobility /x*2 is larger than the self- 
mobility /iii". Furthermore, by considering the special case 




(a) 





(b) 




FIG. 4: Transversal eigenmodes of two trapped spheres (white: re- 
laxed equilibrium state). There are two modes with collective trans- 
lation and relative rotation (a) and two modes with relative translation 
and collective rotation (b). For each mode shown, there is a comple- 
mentary mode with opposite direction of rotation. 



vxx = — M12 1 (for the geometry, see Fig. [2), we infer that 
fi^2~ is positive. Thus, fi^ 1 - > and fi_ < 0, and we find 
A x / 2 > and A 3/i < 0. 

We note that all eigenvalues of mk are positive. To prove 
this statement, we require A„ > in Eqs. ( 127 at and J28at and 
find the condition fi^fi™ 1 - > (/i*^- 1 ) 2 , which is the same 
guaranteeing that the mobility matrix is positive definite. 

The eigenmodes described by Eqs. d27bl i and d28bl > are il- 
lustrated in Fig.|4] The mode in Fig.|4ja) corresponds to the 
eigenvector ai, and the one in Fig.|4jb) to a 3 . Each of the 
modes is complemented by a mode with opposite direction 
of rotation corresponding to the eigenvectors a2 and a.4, re- 
spectively. Fig. |4] shows the faster modes, i.e., Ai > A2 and 
A3 > A4. The reason for this is simply that the relaxations de- 
scribed by ai and a3 show qualitatively the same behavior as 
if the particles were not trapped with respect to their orienta- 
tions (fc r = 0). Hence, a2 and a4 correspond to modes where 
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FIG. 5: Eigenvalues of the transversal modes [Eqs. <27a> and I28al 1 
as a function of the dimensionless center-to-center distance p — 
r/a. The mobilities were calculated using the numerical library HY- 
DROLIB 1 28]. The trap constants are chosen as fc* = fc r . 
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the orientational relaxation is opposite to this "natural" direc- 
tion of rotation. Therefore, these modes are slower compared 
to their respective partners ai and a3. 

These qualitative considerations are in agreement with ex- 
act numbers for the eigenvalues which we plot in Fig. [5] as a 
function of the particle distance. Furthermore, for large dis- 
tances, rotational and transversal motions decouple since A1/3 
tend towards the single-particle translational relaxation rate 
fcV* whereas A 2 /4 assume the corresponding rotational value 
fc r /i r . As obvious from the figure, this decoupling occurs for 
the rotational motion at shorter distances compared to transla- 
tions since rotational flow fields decay faster than translational 
ones. 

IV. TIME-CORRELATIONS IN BROWNIAN MOTION 



J2 n a n <8> b„ = 1. Representing the matrix exponential by 
its spectral decomposition, e~ mkt = ^ n e _A "'a„ (g) b„, we 
obtain from Eq. d32t 

(q m (t + t) q n {t)) = e ~ XlTa i™bin , (33) 

""Tin ^ 

where a/ m (&;„) is the mth (nth) component of the vector 
a.i (hi) and k nn the n-th element of the diagonal matrix k 
(k nn — k 1 , fc', k T , k 1 '). Normalizing the correlation func- 
tion to the square roots of the mean square displacements 
(q^) = k B T jk nn , we finally arrive at 



(q m (t + T)q n (t)) 




To treat the Brownian motion of the two spheres, time- 
dependent random forces and torques f (t) mimicking the mi- 
croscopic degrees of freedom of the surrounding fluid have 
to be added. The total generalized force vector is then f = 
— kq + f (t). Extending Eq. (124b . we obtain the Langevin-type 
equation 

q + mkq = mf (t) , (29) 

which describes a so-called Ornstein-Uhlenbeck process 1 32 ] . 
The random force is assumed to be a Gaussian white noise 
and is fully characterized by its first and second moments, 

(f(t))=0, (f(i) <8>f(t')) = 2k B Tm- 1 S(t-t'). (30) 

The fluctuation-dissipation theorem ( 13 0> relates the second 
moment of the fluctuating forces and torques to the friction 
matrix m _1 13311 . 



The formal solution of the Ornstein-Uhlenbeck process de- 
scribed by Eq. g9j reads I3I 



q(*) = 



-mkt 



q(o) 



dt'e-^-^mf^'), (31) 



where the matrix exponential is defined, as usual, by its Tay- 
lor expansion e _mkt = X]^=o(~ m '^) S / s ' Since we are only 
interested in the fluctuating part of q(t), we omit the deter- 
ministic relaxation e _mkt q(0) due to an initial displacement 
by choosing q(0) = 0. In calculating the correlation matrix, 
we employ the second part of Eq. (I30i and consider the long- 
time limit (t 3> A^ 1 ), where the system has relaxed to thermal 
equilibrium. We finally obtain 



(q(f + t) ® q(f)) = k B Te 



mkri 



r > 



(32) 



Note that the correlation matrix is symmetric, i.e., (q(ii) <8> 

q(«2)) = (q(t 2 )®q(ti)). 

To express the correlation functions explicitly in coordi- 
nates, we introduce the dual eigenvectors b„ via (mk) T b„ = 
A n b ra . Together with the eigenvectors a„, they fulfill the or- 
thonormality and completeness relations, a m • b„ = S mn and 



A. Longitudinal motions 

Longitudinal motions are characterized by the eigenvectors 
and eigenvalues in Eq. i26\ . Translations and rotations are 
decoupled. So the respective auto- (+) and cross- (— ) corre- 
lations for rotational motions are given by [cxp(— /c r /i"'"r) ± 
exp(— /c r /i"'"T)]/2, where = //ji ±A*i2 ■ For translations, 
the correlation functions are the same, however the relaxation 
rates are replaced by fcV±" GHHEJl- 

The correlations for translational motions have already 
been measured experimentally and compared with theoretical 
predictions 1 1 8ll 1 9ll2(JII . These works reveal a good agreement 
between theory and experiment. The correlation functions for 
rotational motions show the same qualitative behavior (Fig.HJl, 
however with some quantitative differences. 

First of all, we realize that the strength of the cross corre- 
lations decreases with increasing particle separation. This is 
explained by the decrease of the cross-mobility /i"^ 1 , which is 




time difference x [arb. units] 

FIG. 6: Correlation functions for longitudinal rotational motions. 
The curve labels indicate the dimensionless center-to-center distance 
p — r/a. The functions are normalized to the mean square angu- 
lar displacement (x| z ) — k^T /k r . The mobilities were calculated 
using the numerical library HYDROLIB 1 28]. 
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due to the spatial decay of the flow fields created by the parti- 
cle rotation. Compared to translational motions, this decrease 
is more pronounced since rotational flow fields decay stronger 
than their translational counterparts. 

Secondly, while the auto-correlation reveals the typical 
monotonous decay, the cross-correlation function features 
an interesting behavior. It is negative (denoted as "anti- 
correlation" in Ref. [18]), since the relative modes decay 
more slowly than the collective ones, as already mentioned 
in Sec. IIII Al Therefore, the negative term in the cross- 
correlation function in the first paragraph dominates. 

The most interesting feature of the cross-correlation func- 
tion is that it exhibits a "memory effect." It vanishes at r = 0, 
in contrast to what one would initially expect for the instanta- 
neous hydrodynamic forces in Stokesian dynamics, and then 
shows a distinct time-delayed extremum (located at r* in 
Fig-EJ- This behavior can be understood as follows. The mo- 
tion of particle 1 creates a fluid flow which instantaneously 
reaches particle 2. However, due to the trap, particle 2 can 
only "react" in a finite time. Thus, the correlation evolves on 
a characteristic time scale which is related to the relaxation 
times and thereby to the trap stiffness. Since the "memory" is 
"stored" in the trap, it cannot last longer than the typical re- 
laxation time. Hence, the correlation decays to zero for times 
larger than r* . 

For sufficiently large particle distances, the self-mobility 
fi"i is much bigger than the cross-mobility /j,"^ ■ Then, to 
leading order, the characteristic time scale is given by the 
single-particle rotational relaxation time r* = (fc r // r ) _1 , in 
analogy to the translational case fl8l Il9l Eoll . Indeed, ac- 
cording to Fig. [6] r* depends only very weakly on the par- 
ticle distance. Nevertheless, the cross-mobility /i 1 ^" is suffi- 
ciently large to separate the two time scales (k T /j, T I ) _1 and 
(k T fi_)~ l , which is the origin of the "anti-correlation". 
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FIG. 7: Correlation functions for transversal motions. The curves 
shown are the mixed cross- and self-correlation for angle and po- 
sition. The curve labels indicate the dimensionless center-to-center 
distance p = r/a. The trap constants are chosen as k = k r . The 
functions are normalized to the square roots of the mean square dis- 
placements (xty) = k B T/k T and (r 2 ix ) = k B T/k t [seeEq. <34ll. 



B. Transversal motions 

The transversal eigenmodes are characterized by Eqs. ( I27> 
and (1281 . Since the two types of eigenvectors are orthogonal 
to each other, i.e., ai/2 • — 0, the dual vectors hx/2 are 
linear combinations of ai and a2 only; the equivalent holds 
for b 3 / 4 . One finds 

^ | a n| 2a n — ( a n ' a ft) a fi (35) 

| a n| 2 | a n| 2 ~~ ( a n ' a ri) 2 

where the index combinations are (n\n) = (1|2), (2|1), (3|4), 
(4 1 3). The correlation functions for the coordinates qi = n x , 
92 = T2x, 93 = Xiy, and 94 = \2y are obtained from Eq. d34l >. 
Now, they are linear combinations of four exponential decays. 

The auto-correlation functions (r\ x (t + T)r\ x {t)) and 
{Xiy(t + T )xiy{t)) show qualitatively the same monotonous 
decay as in the case of longitudinal fluctuations (see Fig.^Jl. 
All the other correlation functions exhibit the features of the 
longitudinal cross-correlation discussed in Sec. 11V Al (except 
for the sign). 



As representative examples, the correlation functions for 
rotations about the y-axis and translations along the x- 
direction, (xiy(t + T)ri x (t)} and (xiy(t + T)r 2x (t)), are plot- 
ted in Fig. In the following, we will refer to them briefly 
as mixed self- and cross-correlation, describing the self- and 
cross-coupling of rotation and translation, respectively. 

Both correlation functions show a distinct time-delayed ex- 
tremum, denoted by r* and r* in Fig. The mixed cross- 
correlation is interpreted in the same fashion as the longitudi- 
nal cross-correlation in the preceeding section. 

The most striking feature is revealed by the mixed self- 
correlation. For a single sphere, translation and rotation are 
not coupled l3lll . i.e., the mixed (self-)correlation vanishes. 
However, in the two-sphere system, there exists a mixed self- 
correlation, as shown in Fig.0 that is mediated by the neigh- 
boring particle. This correlation is weaker than the mixed 
cross-correlation since the flow field created by particle 1 has 
to be reflected by particle 2. As discussed before, particle 2 
"reacts" with a finite delay due to the trap stiffness. Then, in 
addition, particle 1 also has a finite "reaction time." Hence, 
we always expect r * > r* . 

We also studied the influence of the trap stiffnesses fc' and 
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FIG. 8: Mixed self- (label s) and cross-correlation (label c) for var- 
ious trap constant ratios fe : k T . The force constant k was varied 
(fc l = 0.01, 1, and 100), while the torque constant was kept fixed 
(k T = 1). The qualitative behavior of the functions (discussion see 
text) is the same for different particle distances (here: p = 4). 



k r on the correlation functions. As illustrated in Fig. [3] the dif- 
ference in the delay times r* and t* decreases with increasing 
ratio fc' /k" . Furthermore, we observe that the correlations are 
strongest for fc* w fc r . 

Note that the mixed self-correlations cannot be treated 
within the Rotne-Prager approximation. They constitute an 
additional effect of higher order, as mentioned in Sec. Ill CI 



V. CONCLUSIONS 

In this paper, we have presented the complete solution 
for the hydrodynamically coupled translational and rotational 
motions of two colloidal spheres that are harmonically trapped 
with respect to both their positions and orientations. 

Based on pure symmetry arguments and without relying on 
explicit values of the mobilities, we have determined all the 12 
collective eigenmodes and qualitatively discussed their relax- 
ation times. Whereas the properties of the longitudinal modes 
are reminiscent to a system with pure translational degrees of 
freedom, the transversal modes exhibit a characteristic cou- 
pling of translation and rotation. 

In a detailed Langevin-type analysis, we have been able to 
derive the full set of correlation functions characterizing the 
Brownian motion of the particles in the optical traps. The 
analysis relies on the eigenvalue problem of nonsymmetric 



matrices. Explicit examples for the correlation functions at 
different particle separations have been calculated based on 
mobilities which we obtained from the numerical library HY- 
DROLIB lEl . 

The longitudinal fluctuations exhibit the features already 
mentioned in Refs. 0113 m, namely a "memory effect" 
in the cross-correlations. The transversal fluctuations are gov- 
erned by the coupling of translations and rotations. As the 
most striking feature, this coupling is also visible in the self- 
correlations which are governed by a second delay time in ad- 
dition to the one observed in the cross-correlations. The self- 
coupling of translation and rotation for one sphere has to be 
mediated by a second sphere. It is not included in the Rotne- 
Prager approximation and therefore introduces an additional 
effect of higher order. 

Correlation functions involving the rotational degrees of 
freedom are weaker compared to pure translational correla- 
tions since the flow field of a single rotating sphere decays 
as 1/r 2 , compared to translating particles where the decay 
is 1/r. However, the strength of the correlations increases 
for decreasing particle separation which is more pronounced 
whenever rotations are involved. Furthermore, at suffici ently 
small distances, lubrication theory becomes important 13411 . 
and the system introduced in this paper may help to check its 
predictions. 

Our work stresses the rotational degree of freedom and its 
influence on hydrodynamic interactions. We hope to stimu- 
late experimental investigations of the correlation functions 
presented in this paper. 

In the Introduction, we have already mentioned the field of 
microrheology. An extension of our work to viscoelastic me- 
dia based on the theoretical approach presented in Ref. ifnil 
seems appealing. The rotational motion of the probe parti- 
cle results in an interesting deformational mode since it intro- 
duces some torsion in the surrounding medium. It therefore 
can yield additional information about the viscoelastic prop- 
erties. 
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